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The half-filled two-dimensional Hubbard model in presence of a uni- 
form and static parallel magnetic field has been studied by means of the 
Composite Operator Method. A fully self-consistent solution, fulfilling all 
the constrains coming from the Pauli principle, has been found. The rele- 
vant features of a metal-insulator transition in presence of a magnetic field 
have been analyzed. The results qualitatively agree with the ones recently 
obtained by means of experimental investigations. 

PACS numbers: 71.10.Fd,71.27.+a,75.10.-b,71.30.+h 

The response of a two-dimensional (2D) electronic system to a parallel 
magnetic field is very intriguing and several anomalous properties have been 
observed. There is a general agreement that the observed behavior is related 
to the spin polarization, but further studies, both theoretical and experi- 
mental, are needed. In this paper we concentrate on the metal-insulator 
transition (MIT) driven by a in-plane magnetic field. Recent experiments 
on Si-MOSFET [f] and GaAs [2] have shown that by increasing field the 
spin system polarizes and the system undergoes a MIT before reaching the 
full polarization. Apparently, an important role is played by the electron- 
electron interaction, being r s = U/K (the ratio of Coulomb interaction 
energy to the mean kinetic energy) very large. 

In order to make a qualitative and preliminary study of this phenomenon 
we consider the 2D Hubbard model in presence of a parallel external mag- 
netic field. Since a parallel field does not couple to the orbital motion of 
electrons, the Hamiltonian is given by 

H = J^(-4ta ii -Li5 ii )cHi) c(j) + C/^n T (i) n 4 (i) - \h £ n 3 (i) (1) 

ij i i 

where c (i) and c* (i) are the annihilation and creation operators of electrons 
in spinorial notation; i = (i, t) where i are vectors of a 2D Bravais lattice; /i 
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is the chemical potential; a^j denotes the projector on first-neighbor sites; 
U is the local Coulomb interaction, n a (i) = c a (i) c a (i) is the charge den- 
sity of the electrons with spin a; n 3 (i) is the third component of the spin 
density operator; h is proportional to the intensity of the external magnetic 
field. In the framework of the Composite Operator Method (COM) [3], we 
introduce the basis ip^(i) = (^(i),vHi)) where = (1 — n(i))c(i) and 
rj(i) = n(i) c(i) are the Hubbard operators responsible for the transitions 
|0) <-> \a) and \a) <-> |||), respectively. The composite operator ip(i) satisfies 
the equation of motion 

ij^ (*) = (\ha 3 - 1> (*)-24 (1 + r 3 ) c a (1 - r 3 ) 77 (i)-4t r 3 ir (i) 

where cj acts on the spin degree of freedom a =|, j and r on the internal 
degree of freedom ip = ^ , rj. a and r are Pauli matrices. We also use the no- 
tation 4> a (i, t) = £\ a-ycf) (j, Moreover, we have 7r(z) = ^cj m c a (i) + 
[c ta (i) f](i)] where cr M = (1,<t), cr M = (-1,<?) and n M (z) = 0^(1)0^0(1) 
describe the total charge- (ji = 0) and spin- (ji = 1, 2, 3) density operators. 

In the polar approximation [3] we linearize the equation of motion by 
projecting the source on the basis ip(i). Then, the retarded Green's function 
<S(k, lo) = T (1Z [ip(i) ip' (j)] ), where T and 1Z are the Fourier transform and 
the usual retarded operators, respectively, has the following expression 

where the energy spectra (k) are the eigenvalues of the energy ma- 
trix e (k) = T ({ J (i, t) , ip^ (j, t) }) / _1 (k) and the spectral density matrices 

erW (k) are calculated by means of the formula <r^(k) = Q a i(\z) ^^(k) I^p 
where f2(k) is the matrix whose columns are the eigenvectors of the energy 
matrix e (k) and /(k) = T ({V> (i, t) , ^ (j, t)}) is the normalization matrix. 
The explicit expressions of (k) and (k) will be given elsewhere. 
Calculations show that the Green's function depends on the following set 
of parameters: /j, m, A CT , p a . m = ^ (713 (i)) is the magnetization per site. 
The parameters A CT and p a describe a constant shift of the bands and a 
band width renormalization, respectively, and are defined as 

A* = (£(0&(0)-(t£(0»£(0) W 
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P. = ^ [<r£ (<) (<)> + 2 (-)* K (<) n 3 (7)}] - (ft (i) ^ r/J (i) ^ (t) 

(5) 




Fig. 1. (left) The magnetization m as a function of the external magnetic field h for 
T = 0, n = 1 and various values of the Coulomb repulsion U ; (right) The critical 
value U c of the Coulomb repulsion U for the MIT as a function of the external 
magnetic field h for T = and n = 1. t is taken as unity. 



The determination of these parameters is very crucial and wrong re- 
sults are easily obtained as shown in Ref. [4]. The parameters m and 

are expressed in terms of the Green's function as m — ^ (C44 — C22)) 
Af = Cf x — C22 and Aj = C33 — C44. We have defined the correlation ma- 
trices C = (ip (i) ip> (z)) and C a = (ip a (i) ifi* (i)). The other parameters /i 
and p a are not determined by the equation of motion and are fixed by choice 
of the representation where the Green's functions are realized [5]. In the 
COM we choose the representation by requiring that all the relations among 
the operators dictated by the algebra (Pauli principle) are conserved also at 
the level of expectation values. In the present study, this requirement leads 
to Cn = C33 and C12 = C34 = 0. Because we are interested in the study 
of the MIT, we consider the special case of half filling (n = {n(i)) = 1) 
where: /i = ^, Aj = — A|, P] = Pi — 2m and C12 = C34 = 0. It is worth 
to note that these latter relations are a manifestation of the particle-hole 
symmetry which is conserved owing to the choice of the representation. 
Any other choice of the representation will lead, in the context of the pole- 
approximation, to a violation of the symmetry [4]. Finally, we have a set of 
three coupled self-consistent equations which determine the three parame- 
ters which are left: m, A = At, p = p-t. 

In Fig. 1 (left panel) we plot the magnetization m versus the magnetic 
field h. The magnetization is an increasing function of both the applied 
magnetic field and the Coulomb interaction U. It reaches the saturation 
value (i.e., 1/2) at a critical value of the magnetic field, which depends on 
the intensity of the Coulomb interaction. At zero temperature T = 0, when 
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U approaches the bandwidth Wzd = 8i, the magnetization experiments a 
discontinuous jump from zero up to the saturation value. We have also 
inspected the analytical behavior of the static susceptibility by analyzing 
the self-consistent equations in the limit of very low magnetic fields. Results 
show a divergence when the Hubbard repulsion approaches the bandwidth 
at zero temperature. The double occupancy decreases when increasing both 
the interaction and the magnetic field. The latter provides the spins of the 
electrons with an orientation and, due to the Pauli principle, reduces the 
double occupancy. There is a quite good agreement between COM results 
and Gutzwiller ones [6]. 

The MIT can be studied by looking at the density of states (DOS): the 
opening of a gap in the DOS is a signal of the transition from metallic to 
insulating phase. In Ref. [7] we have studied the MIT exhibited by the Hub- 
bard model in absence of magnetic field for the 2D and 3D cases. It was 
found that the transition is driven by the Coulomb interaction: there is a 
critical value U c where the MIT occurs. In particular, the value U c = 1.68W 
(W2D = 8t and Wzd = 12i for the 2D and 3D system, respectively) was 
reported. In presence of a magnetic field the value of U c is drastically influ- 
enced. In Fig. 1 (right panel) we plot the critical value versus the magnetic 
field at zero temperature. As we turn on a rather small magnetic field, 
the critical value U c suddenly jumps from U c = 1.68W to U c = W. This 
discontinuity at zero field is related to the discontinuity of the magnetiza- 
tion, as shown in Fig. 1 (left panel). By increasing h, U c decreases and 
vanishes when the field equates the bandwidth at zero temperature (i.e., 
U c (h,T = 0) = W — h), in qualitative agreement with the experimental 
findings. 

In conclusion, our study shows that the 2D Hubbard model in presence 
of a parallel magnetic field can describe the experimental evidence of a field- 
driven MIT. The transition is controlled by the field and disappears for some 
critical value of it. A more detailed discussion of the MIT and of the order 
parameter controlling the transition will be reported elsewhere. 
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